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[0001] REDUCED COMPLEXITY SLIDING WINDOW BASED EQUALIZER 

[0002] CROSS REFERENCE TO RELATED APPLICATION(S) 

[0003] This application claims priority from U.S. provisional application no. 

60/452,165, filed on March 3, 2003, which is incorporated by reference as if fully set 

forth. 

[0004] FIELD OF INVENTION 

[0005] The invention generally relates to wireless communication systems, In 
particular, the invention relates to data detection in such systems. 

[0006] BACKGROUND 

[0007] Due to the increased demands for improved receiver performance, many 

advanced receivers use zero forcing (ZF) block linear equalizers and minimum mean 
square error (MMSE) equalizers. 

[0008] In both these approaches, the received signal is typically modeled per 

Equation 1. 

r = Hd + n 

Equation 1 

[0009] r is the received vector, comprising samples of the received signal. H is 

the channel response matrix, d is the data vector. In spread spectrum systems, such 
as code division multiple access (CDMA) systems, d is the spread data vector. In 
CDMA systems, data for each individual code is produced by despreading the estimated 
data vector d with that code, n is the noise vector. 

[0010] In a ZF block linear equalizer, the data vector is estimated, such as per 

Equation 2. 

d = (H)'r 
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Equation 2 

[0011] Q H is the complex conjugate transpose (or Hermetian) operation. In a 

MMSE block linear equalizer, the data vector is estimated, such as per Equation 3. 
d = (H H H + a 2 l) 1 r 

Equation 3 

[0012] In wireless channels experiencing multipath propagation, to accurately 

detect the data using these approaches requires that an infinite number of received 
samples be used. One approach to reduce the complexity is a sliding window approach. 
In the sliding window approach, a predetermined window of received samples and 
channel responses are used in the data detection. After the initial detection, the 
window is slid down to a next window of samples. This process continues until the 
communication ceases. 

[0013] By not using an infinite number of samples, an error is introduced into 

the data detection. The error is most prominent at the beginning and end of the 
window, where the effectively truncated portions of the infinite sequence have the 
largest impact. One approach to reduce these errors is to use a large window size and 
truncate the results at the beginning and the end of the window. The truncated 
portions of the window are determined in previous and subsequent windows. This 
approach has considerable complexity. The large window size leads to large 
dimensions on the matrices and vectors used in the data estimation. Additionally, this 
approach is not computationally efficient by detection data at the beginning and at the 
ends of the window and then discarding that data. 

[0014] Accordingly, it is desirable to have alternate approaches to data detection. 

[0015] SUMMARY 

[0016] Data estimation is performed in a wireless communications system. A 

received vector is produced. For use in estimating a desired portion of data of the 
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received vector, a past, a center and a future portion of a channel estimate matrix is 
determined. The past portion is associated with a portion of the received signal prior to 
the desired portion of the data. The future portion is associated with a portion of the 
received vector after the desired portion of the data and the center portion is associated 
with a portion of the received vector associated with the desired data portion. The 
desired portion of the data is estimated without effectively truncating detected data. 
The estimating the desired portion of the data uses a minimum mean square error 
algorithm having inputs of the center portion of the channel estimate matrix and a 
portion of the received vector. The past and future portions of the channel estimate 
matrix are used to adjust factors in the minimum mean square error algorithm. 

[0017] BRIEF DESCRIPTION OF THE DRAWING(S) 

[0018] Figure 1 is an illustration of a banded channel response matrix. 

[0019] Figure 2 is an illustration of a center portion of the banded channel 

response matrix. 

[0020] Figure 3 is an illustration of a data vector window with one possible 

partitioning. 

[0021] Figure 4 is an illustration of a partitioned signal model. 

[0022] Figure 5 is a flow diagram of sliding window data detection using a past 

correction factor. 

[0023] Figure 6 is a receiver using sliding window data detection using a past 
correction factor. 

[0024] Figure 7 is a flow diagram of sliding window data detection using a noise 

auto-correlation correction factor. 

[0025] Figure 8 is a receiver using sliding window data detection using a noise 

auto-correlation correction factor. 
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[0026] DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT(S) 
[0027] Hereafter, a wireless transmit/receive unit (WTRU) includes but is not 

limited to a user equipment, mobile station, fixed or mobile subscriber unit, pager, or 
any other type of device capable of operating in a wireless environment. When referred 
to hereafter, a base station includes but is not limited to a Node-B, site controller, 
access point or any other type of interfacing device in a wireless environment. 
[0028] Although reduced complexity sliding window equalizer is described in 
conjunction with a preferred wireless code division multiple access communication 
system, such as CDMA2000 and universal mobile terrestrial system (UMTS) frequency 
division duplex (FDD), time division duplex (TDD) modes and time division 
synchronous CDMA (TD-SCDMA), it can be applied to various communication system 
and, in particular, various wireless communication systems. In a wireless 
communication system, it can be applied to transmissions received by a WTRU from a 
base station, received by a base station from one or multiple WTRUs or received by one 
WTRU from another WTRU, such as in an ad hoc mode of operation. 
[0029] The following describes the implementation of a reduced complexity 

sliding window based equalizer using a preferred MMSE algorithm. However, other 
algorithms can be used, such as a zero forcing algorithm. h( ) is the impulse response 
of a channel. d{k) is the k th transmitted sample that is generated by spreading a 

symbol using a spreading code. It can also be sum of the chips that are generated by 
spreading a set of symbols using a set of codes , such as orthogonal codes, r(-) is the 
received signal. The model of the system can expressed as per Equation 4. 

oo 

r(r)= ^d(k)h(t-kT c )^n(t) -oo<r<oo 

Equation 4 

[0030] n(t) is the sum of additive noise and interference (intra-cell and inter-cell). 

For simplicity, the following is described assuming chip rate sampling is used at the 
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receiver, although other sampling rates may be used, such as a multiple of the chip 
rate. The sampled received signal can be expressed as per Equation 5. 



r(j)= ^d(k)h(j-k) + nU)= ^dU-k)h(k) + n(j) je {...,-2,-1,0,1,2,...} 

*=-~ *=— 

Equation 5 

T c is being dropped for simplicity in the notations. 

[003 1] Assuming h( ) has a finite support and is time invariant. This means that 

in the discrete-time domain, index L exists such that h(i) = 0 for i < 0 and i > L . As a 
result, Equation 5 can be re-written as Equation 6. 

r(j) = ^h(k)d(j-k) + n(j) je {...,-2,-1,0,1,2,...} 

Equation 6 

[0032] Considering that the received signal has M received signals 

r(0), • • • , r(M - 1) , Equation 7 results. 
r = Hd + n 

where, 

r = [r(0),-,r(M-l)] r eC M , 

d = [d(-L + l),d(-L + 2),...,d(0),d(l),...,d(M -1)] T e C M+L=l 
n=[n(0),"-,n(M-l)] 7 "e C M 



H = 



h(L-l) h(L-2) 

0 h(L-l) h(L-2) 



0 



h(l) HO) 0 

h(l) h(0) 0 



e £Mx(M+L-l) 



h(L-l) h(L-2) ■■■ h(i) /i(0)_ 

Equation 7 

[0033] Part of the vector d can be determined using an approximate equation. 

Assuming M > L and defining N = M - L + 1, vector d is per Equation 8. 

d = [d(-L + l),tf(-L + 2),...,rf(-l),rf(0),rf(l),...,rf(N - l) y d(N) y ...,d(N + L - 2)] T e C N+2L ~ 2 



L-\ 



L-l 
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Equation 8 

[0034] The H matrix in Equation 7 is a banded matrix, which can be represented 

as the diagram in Figure 1. In Figure 1, each row in the shaded area represents the 
vector [/i(L-l),ft(L-2),...,/i(l),ft(0)], as shown in Equation 7. 

[0035] Instead of estimating all of the elements in d, only the middle N elements 

of d are estimated, d is the middle N elements as per Equation 9. 

a = [d(0) y ...,d(N-i)] T 

Equation 9 

[0036] Using the same observation for r, an approximate linear relation between 

r and d is per Equation 10. 
r = Hd + n 

Equation 10 

[0037] Matrix H can be represented as the diagram in Figure 2 or as per 

Equation 11. 



H = 



h(0) 


0 




h(l) 


h(0) '■ 






h(l) ■■ 


0 


h(L-l) 




HQ) 


0 


h(L-l) '■ 


h(l) 



0 



h(L-l) 



Equation 11 

[0038] As shown, the first L-l and the last L-l elements of r are not equal to the 

right hand side of the Equation 10. As a result, the elements at the two ends of vector 
d will be estimated less accurately than those near the center. Due to this property, a 
sliding window approach is preferably used for estimation of transmitted samples, such 
as chips. 

[0039] In each, k th step of the sliding window approach, a certain number of the 

received samples are kept in r [k] with dimension iV+L-1. They are used to estimate a 
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set of transmitted data d[k] with dimension N using equation 10. After vector d[k] is 

estimated, only the "middle" part of the estimated vector d[k] is used for the further 

data processing, such as by despreading. The "lower" part (or the later in-time part) of 

d[k] is estimated again in the next step of the sliding window process in which r [k+1] 

has some of the elements r [k] and some new received samples, i.e. it is a shift (slide) 
version of r [k] . 

[0040] Although, preferably, the window size N and the sliding step size are 

design parameters, (based on delay spread of the channel (L), the accuracy requirement 
for the data estimation and the complexity limitation for implementation), the 
following using the window size of Equation 12 for illustrative purposes. 
N = 4N s x SF 

Equation 12 

SF is the spreading factor. Typical window sizes are 5 to 20 times larger than the 
channel impulse response, although other sizes may be used. 

[0041] The sliding step size based on the window size of Equation 12 is, 

preferably, 2N s xSF. e {1,2,...} is, preferably, left as a design parameter. In 
addition, in each sliding step, the estimated chips that are sent to the despreader are 
2N S x SF elements in the middle of the estimated d[k] . This procedure is illustrated in 
Figure 3. 

[0042] One algorithm of data detection uses an MMSE algorithm with model 

error correction uses a sliding window based approach and the system model of 
Equation 10. 

[0043] Due to the approximation, the estimation of the data, such as chips, has 

error, especially, at the two ends of the data vector in each sliding step (the beginning 
and end). To correct this error, the H matrix in Equation 7 is partitioned into a block 
row matrix, as per Equation 13, (step 50). 

h = [h, |h |hJ 
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Equation 13 

[0044] Subscript "p" stands for "past", and "/" stands for "future". H is as per 

Equation 10. H p is per Equation 14. 



h(L-l) h(L-2) 

0 h(L-l) 

0 
0 

0 



- h(l) 

- h(2) 

0 h(L-l) 
0 

0 



eC 



(N+L-i)x(L-l) 



[0045] 



H 7 is per Equation 15. 

0 0 



(Af+£,-l)X(L-l) 



Equation 14 



0 0 

h(0) 0 ••• 0 eC 

! •. '•. 0 

h(L-3) ■■■ h(0) 0 

h(L-2) h(L-3) ••• /i(0) 

Equation 15 

[0046] The vector d is also partitioned into blocks as per Equation 16. 

IT l AT 



d=k,|a r Kr 



Equation 16 



[0047] d is the same as per Equation 8 and d p is per Equation 17. 

d =[d(-L+l) d(-L + 2) ••• d(-l)] T 6C t_1 



Equation 17 



[0048] d f is per Equation 18. 
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d f =[d(N) d(N + l) - d(N + L-2)J e C L ~ l 

Equation 18 

[0049] The original system model is then per Equation 19 and is illustrated in 

Figure 4. 

r = H p d p +Hd + H / d / +n 

Equation 19 

[0050] One approach to model Equation 19 is per Equation 20. 

r = Hd + n x 

where r =r-H p d p and n 1 = H / d / + n 

Equation 20 

[0051] Using an MMSE algorithm, the estimated data vector d is per Equation 
21. 

a = g rf H w (^HH w +2 1 )- , F 

Equation 21 

[0052] In Equation 21, g d is chip energy per Equation 22. 

E{d(i)d\j)}= gAj 

Equation 22 

[0053] r is per Equation 23. 

Equation 23 

[0054] d p , is part of the estimation of d in the previous sliding window step. J L l 

is the autocorrelation matrix of n, , i.e., Z t = Ejnjn^}. If assuming H / d / and n are 

uncorrelated, Equation 24 results. 
^= 8d U f H H f ^E{nn H } 

Equation 24 
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[0055] The reliability of d p depends on the sliding window size (relative to the 

channel delay span L) and sliding step size. 

[0056] This approach is also described in conjunction with the flow diagram of 

Figure 5 and preferred receiver components of Figure 6, which can be implemented in a 
WTRU or base station. The circuit of Figure 6 can be implemented on a single 
integrated circuit (IC), such as an application specific integrated circuit (ASIC), on 
multiple IC's, as discrete components or as a combination of IC('s) and discrete 
components. 

[0057] A channel estimation device 20 processes the received vector r producing 

the channel estimate matrix portions, H p , H and U f , (step 50). A future noise auto- 
correlation device 24 determines a future noise auto-correlation factor, g d H f H" , (step 

52). A noise auto-correlation device 22 determines a noise auto-correlation factor, 
isjnn" }, (step 54). A summer 26 sums the two factors together to produce Z l , (step 56). 
[0058] A past input correction device 28 takes the past portion of the channel 

response matrix, H p , and a past determined portion of the data vector, d , to produce 

a past correction factor, H p d p , (step 58). A subtracter 30 subtracts the past correction 

factor from the received vector producing a modified received vector, r , (step 60). An 
MMSE device 34 uses Sj , H , and r to determine the received data vector center 

A 

portion d , such as per Equation 21, (step 62). The next window is determined in the 

A 

same manner using a portion of d as d p in the next window determination, (step 64). 

As illustrated in this approach, only data for the portion of interest, d , is determined 
reducing the complexity involved in the data detection and the truncating of unwanted 
portions of the data vector. 

[0059] In another approach to data detection, only the noise term is corrected. In 

this approach, the system model is per Equation 25. 
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r = Hd + n 2 , where n 2 = H p d p +H / d / +n 

Equation 25 

[0060] Using an MMSE algorithm, the estimated data vector d is per Equation 

26. 

S = g (i H ff (^HH ,, +S 2 )- , r 

Equation 26 

[0061] Assuming H p d p , H / d / and n are uncorrected, Equation 27 results. 
X 2 =g d K p n H p +g d H f ll H f +E{nn H } 

Equation 27 

[0062] To reduce the complexity in solving Equation 26 using Equation 27, a full 

matrix multiplication for H p H" and I^H" are not necessary, since only the upper 

and lower corner of H p and H f , respectively, are non-zero, in general. 

[0063] This approach is also described in conjunction with the flow diagram of 

Figure 7 and preferred receiver components of Figure 8, which can be implemented in a 
WTRU or base station. The circuit of Figure 8 can be implemented on a single 
integrated circuit (IC), such as an application specific integrated circuit (ASIC), on 
multiple IC's, as discrete components or as a combination of ICCs) and discrete 
components. 

[0064] A channel estimation device 36 processes the received vector producing 
the channel estimate matrix portions, H^HandH,, (step 70). A noise auto- 
correlation correction device 38 determines a noise auto-correlation correction factor, 
g d H p H H p + g rf H r H " , using the future and past portions of the channel response matrix, 

(step 72). A noise auto correlation device 40 determines a noise auto-correlation factor, 
£{nn"}, (step 74). A summer 42 adds the noise auto-correlation correction factor to the 
noise auto-correlation factor to produce S 2 , (step 76). An MMSE device 44 uses the 
center portion or the channel response matrix, H, the received vector, r , and Z 2 to 
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estimate the center portion of the data vector, d , (step 78). One advantage to this 
approach is that a feedback loop using the detected data is not required. As a result, 
the different slided window version can be determined in parallel and not sequentially. 
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